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Abstract 

Observational data hint at a finite universe, with spherical manifolds 
such as the Poincare dodecahedral space tentatively providing the best 
fit. Simulating the physics of a model universe requires knowing the 
eigenmodes of the Laplace operator on the space. The present article 
provides explicit polynomial eigenmodes for all globally homogeneous 
3-manifolds: the Poincare dodecahedral space S 3 /I*, the binary oc- 
tahedral space S 3 /0*, the binary tetrahedral space S 3 /T*, the prism 
manifolds S 3 /D^ n and the lens spaces L(p, 1). 

1 Introduction 

The past decade has seen intense work on multiconnected 3-manifolds as 
models for the physical universe. Well-proportioned 3-manifolds explain the 
missing broad-scale fluctuations in the cosmic microwave background, first 
discovered by the COBE satellite [I and later confirmed by the WMAP 
satellite [21 El- (Please see Ref. jlj for an elementary exposition.) A well- 
proportioned 3-manifold is one whose three dimensions are of comparable 
magnitudes. Ill-proportioned manifolds, with one dimension significantly 
larger or smaller than the other two, fail to explain the missing broad-scale 
fluctuations and indeed predict exactly the opposite, namely elevated broad- 
scale fluctuations [S], contrary to observations. Current work, therefore, 
focuses on well-proportioned spaces. 

The density of ordinary matter alone would suggest a hyperbolic universe, 
and so ten years ago researchers studied hyperbolic models. The situation 
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changed dramatically in late 1998 with the discovery of a still-mysterious 
vacuum energy that raises the universe's mass-energy density parameter Q to 
the level required for a flat space (Q = 1) or a slightly spherical space (Q > 1). 
The first-year WMAP paper put Q at 1.02 ±0.02 at the la level, while the 
three-year WMAP paper [3] reports six different distributions ranging from a 
rather flat fl = 1.0031q;^3 ^° a surprisingly curved Q = 1.037^q;q2i depending 
on what external data set one uses to resolve the geometrical degeneracy. 

The Poincare dodecahedral space [B], defined as the quotient S 3 /I* of 
the 3-sphere by the binary icosahedral group I*, explains both the missing 
fluctuations and the observed mass-energy density |7] and so researchers are 
now modelling it more precisely for better comparison to observations. State- 
of-the-art simulations find that the Poincare dodecahedral space matches 
observed broad-scale fluctuations when 1.015 < Q < 1.020 jH El or 1.022 < 
Q < 1.034 |10j . in excellent agreement with the three-year WMAP estimates 
of fl However, other topologies, notably the quotient S 3 /0* of the 3-sphere 
modulo the binary octahedral group, also remain viable [9*1 110j. 

Attempts to confirm a multiconnected universe using the circles-in-the- 
sky method [TT] have failed [12] ■ (Exception: one research group claims to 
have found hints of matching circles [13] and a second group has confirmed 
the match [T3j, but the former provide no statistical analysis while the lat- 
ter find the match to lie below the false positive threshold, so the former's 
claim remains unconvincing.) Accepting the result that any potential circle 
pairs are undetectable [T2J, the question remains: are the circles really not 
there, or are they merely hidden by various sources of contamination, such as 
the Doppler and integrated Sachs- Wolfe components of the microwave back- 
ground? Answering this question requires great care, because the level of 
contamination depends both on the topology and on the choice of cosmolog- 
ical parameters. Results so far remain inconsistent: the circle-searchers' own 
analysis finds their negative results to be robust for a dodecahedral universe 
in spite of the contamination [Tl], while other researchers find contamination 
strong enough to hide matching circles [To"| . 

To determine the observational consequences of a given cosmological 
model, researchers take a Fourier approach and express physical quantities, 
such as density fluctuations in the primordial plasma, as linear combinations 
of the eigenmodes of the Laplacian (more briefly, the modes), which can then 
be integrated forward in time. Thus all studies of cosmic topology require 
knowing the modes of the 3-manifold under consideration. Different research 
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groups have determined the modes in different ways (HI El E3 EH HZ| , but so 
far all approaches have required some sort of numerical computation, either 
the extraction of the eigenvectors of a large matrix or the solution of a large 
set of simultaneous equations. 

The present article provides the modes as explicit polynomials with in- 
teger coefficients. That is, for each manifold S 3 /T and each wavenumber k, 
we provide a finite set of T-invariant polynomials of degree k spanning the 
full space of modes. We provide these polynomials for the binary icosahedral 
space S 3 / 1* (better known as the Poincare dodecahedral space), the binary 
octahedral space S 3 /0*, the binary tetrahedral space S 3 /T*, the binary di- 
hedral spaces S 3 /D^ n (better known as prism manifolds), the homogeneous 
lens spaces L(p, 1), and the 3-sphere itself. These spaces comprise the full set 
of globally homogeneous spherical 3-manifolds (called single action spaces in 
the classification of Ref. [TH]). 

The ideas in the present article draw heavily on Klein's 1884 Vorlesun- 
gen iiber das Ikosaeder JUL extending Klein's work to produce full bases of 
explicit polynomials. 

Sections 121 through El are elementary in nature, laying a foundation, es- 
tablishing terminology, and translating into geometric language some con- 
cepts that have recently appeared in the cosmic topology literature only in 
quantum mechanical bracket language. Sections through provide the 
real content of this article, namely the explicit polynomials for the modes 
of S 3 /I\ S 3 /0\ S 3 /T\ S 3 /D* m , L(p,l) and S 3 . Section M outlines future 
work. 

2 Lifting S 2 to S 3 

Let us parameterize the 2-sphere S 2 first as the Riemann sphere C U {oo} 
and then, more elegantly, as the complex projective line CP 1 . First map the 
unit 2-sphere onto the equatorial plane via stereographic projection from the 
south pole: typical points map via (x,y,z) i— ► + i-^ while the south 
pole (0,0,-1) maps to oo. To accommodate the south pole into the same 
format as the typical points, write each image point as a formal fraction. 
Each typical point maps to a formal fraction while the south pole maps 
to the formal fraction |. Two formal fractions § and are equivalent if 
and only if a[3' = a' (3. Let [%] denote the equivalence class of all formal 
fractions ^ equivalent to §. The set of all such equivalence classes defines 
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Figure 1: The dihedron, tetrahedron, octahedron and icosahedron project 
radially to tilings of the 2-sphere. The dihedron has only two faces, one 
filling the whole northern hemisphere and the other filling the whole southern 
hemisphere. The dihedron's vertices and edges lie evenly along the equator. 
(Artwork by Adam Weeks Marano.) 

the complex projective line CP 1 . The indeterminate fraction jj is of course 
excluded from the discussion. In summary, the set CP 1 of equivalence classes 
[^] of of nontrivial complex formal fractions parameterizes the 2-sphere S 2 . 

Writing the formal fraction ^ as an ordered pair (a, (3) G C 2 immediately 
yields a map C 2 — (0,0) — > CP 1 sending each point (a,/3) G C 2 to the 
corresponding equivalence class [§] G CP 1 . The radial direction in C 2 is 
largely irrelevant, because (a,/?) and (ra,r/?) map to the same class [§] = 
G CP 1 for all real r ^ 0. We may therefore restrict the map's domain to 
the unit 3-sphere S 3 C C 2 without compromising its image in CP 1 . Recalling 
from the preceding paragraph that CP 1 parameterizes S 2 , this yields a map 
S 3 ^S 2 . 

The Clifford flow on S 3 G C 2 is the fixed-point-free motion taking (a, (3) \— > 
(e lt a,e lt f5) for < t < 2n. Under the Clifford flow, each point (a, (3) traces 
out a great circle, which collectively form a set of Clifford parallels. The 
points along a given Clifford parallel all map to the same equivalence class 
[f^f ] = [f] e CP 1 . In other words, the preimage of each point of S 2 is a 
Clifford parallel in S 3 which we will call a fiber of the map. 

3 Constructing the Groups 

The classic platonic solids project radially to tilings of S 2 (Figure HJ). We 
ignore the cube and dodecahedron, because they are dual to the octahedron 
and icosahedron respectively, but include the dihedron, which works out fine 
as a tiling of S 2 even though it falls flat as a traditional polyhedron. The 
orientation-preserving symmetries of these tilings comprise the dihedral group 
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D m , the tetrahedral group T, the octahedral group O and the icosahedral group 
I, of order 2m, 12, 24 and 60, respectively. 

Every matrix ( a c h d ) G GL 2 C respects Clifford parallels in the sense 
that left-multiplication takes each fiber (e**a, e lt f3) to another fiber {e lt {aa + 
bf3),e lt (ca + d[3)). Therefore the matrix action projects down to a well- 
defined map S 2 — > S 2 taking p|] to [ ^^f ]■ Restricting our attention to 
unitary matrices ( ° d ) G U(2) ensures rigid motions of C 2 , which in turn 
project down to rigid motions of S 2 . The unitary matrices, however, still al- 
low room for "sliding along the fibers", for example via the matrix ( q ° eie ), 
which of course has no effect when projected down to S 2 . To obtain an (al- 
most) unique matrix for each rotation of S 2 , restrict to the special unitary 
group SU(2), whose matrices take the form ( * ~™ ). Each rotation of S 2 is 
realized by exactly two special unitary matrices ±( z w _ f ). 

How may we construct the pair of matrices ±( z w ~™ ) realizing a given 
rotation of S 2 7 To construct an order n rotation about a desired fixed point 

we require that 

The phase factor e~ ensures that we rotate the correct amount about the 
given fixed point, so the isometry will have the desired order n. A quick 
calculation gives 

f z \ = ± 1 / p q \ ( e^p 
\ w J pp + qq\-qp)y e^rq 

and ±( ^ ~™ ) follows immediately. 

Equation (j2J) makes it easy to write down matrices for the groups D m , 
T, O and /. Two matrices realize each rotation of S 2 , so for example the 
4-element group D 2 C SO(3) is realized by an 8-element matrix group D\ C 
SU(2) called the binary dihedral group of order 8 (Table Q). One might 
hope to extract a 4-element subgroup of D\ realizing D2 directly, but this 
is impossible because squaring an element of order 2 always gives ( _ o -1 )> 
never (J5), and once ( _ ) is in a subgroup, so is the negative of every 
matrix in that subgroup. 

Still using Equation @, one may easily write down matrices for the 12- 
element tetrahedral group T C 50(3), giving the 24-element binary tetrahe- 
dral group T* C SU(2), and similarly for the 48-element binary octahedral 
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center 

0, oo 
±1 
±i 



order 
1 

2 

2 

2 



matrices 




Table 1: The four rotations in D2 C SO (3) correspond to eight matrices in 
D* C SU(2). 

group O* and the 120-element binary icosahedral group I*. The quotient 
S 3 /I* defines the Poincare dodecahedral space. 



4 Constructing Symmetric Polynomials 

Parameterize the 2-sphere S 2 as CP 1 (Section EJ and let P = {[g-], . . . , [£a]} 
be a set of points thereon, whose symmetry we hope to capture in a polyno- 
mial. As a starting point, the polynomial 

(*-*)... (3) 

qi q n 

has roots exactly at P. Replacing the variable z with a formal fraction 
[|] e CP 1 gives 

(|-^)-(^-^) (4) 
P qi P q n 



which Klein, in his 1884 Vorlesungen iiber das Ikosaeder [19] . writes as a 
homogeneous polynomial 

(qia-pxP)---(q n a-p n /3). (5) 



I thank Peter Kramer for his recent article [20] pointing out the relevance of 
Klein's work to current investigations in cosmic topology. 
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To test whether the polynomial JSJ) is invariant under a symmetry 7 = 
( z m ~™ ) G 5'f/(2), rewrite expression © as 

n *)(; -;)(;) c) 

and let 7" 1 = ( ™ ) act on (a, (3), transforming the polynomial to 

n 

n 

i=l 

which equals 

n 

8=1 

In effect the action of 7 transforms each root (pi, qj) according to the rule 

(ft,ft)->(ft,?i)(4y) = (7)(s)- 

Example SJl. Consider the points P = {y,y,— j-,^-} with polynomial 

(a -p){pt- i/3)(a +P)(a + i(3) = a 4 - f3\ (9) 

and consider the matrix 72 = ( q _° ) G SU(2). Geometrically, 72 projects 
down to an order 2 rotation of S 2 about the north pole, which interchanges 
the points of P in pairs. Acting on S 3 C C 2 , 72 takes the roots to 

r i — 1 — i 1, — 1 — i 1 i, . , 

i p = <=?• -v -} = {-'-■ -7- < 10 > 

with formal fractions treated as vectors, so for example — i— j = . The 

action of 72 permutes the roots and also multiplies them by a factor of —i. 
Fortunately these changes leave the polynomial invariant, 

(-i)(a + p) (-i)(a + z/3) (-i)(a - /3) (-i)(a - = a 4 - [3 4 . (11) 

Example |4j2. Consider the same points P = {j, y, -y, — y} with the same 
polynomial a 4 — /3 4 as in Example EJl, but now act by 74 = ( e ^ _° <7r/4 ). 



( Pi % ) 



2; w 
-w z 



a 
P 



(7) 



( Pi Qi) 



z w 
-w z 



(8) 
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Geometrically, 74 projects down to an order 4 rotation of S 2 about the north 
pole, which permutes the points of P cyclically. Acting on S 3 C C 2 , 74 takes 
the roots to 

74 P = {e~ i7v/4 j, e~ in/4 -^, e^ /4 ^-, e" m/4 i}. (12) 

Again the roots have been permuted, but this time they are multiplied by a 
factor of e _Mr//4 , so the polynomial maps to 

(e" i7r/4 ) 4 (a - ip){a + p){a + i[3){a -(3) = -(a 4 - /? 4 ). (13) 

Thus 74 does not leave this polynomial invariant, but rather sends it to its 
negative. 



5 Generalized Complex Derivatives 

Let us extend the complex derivative operator from the class of complex- 
differentiable (i.e. analytic) functions to the broader class of real-differentiable 
(i.e. smooth) functions. 

A complex- valued function / of a complex variable z = x + iy is differen- 
tiable in the complex sense if and only if it is different iable in the real sense 
( as a function of x and y) and, in addition, 

dy % dx' ^ 

If Equation ()14|) holds, the complex derivative || is defined to be the common 
value of |^ and 4|^. If Equation does not hold, then the function / 
has no complex derivative. For example, complex conjugation f(z) — z fails 
test (|14|) and the symbol |j has no meaning in this context. 
An alternative definition of the derivative operator 

8 1/9 +i^ (15) 



dz 2 \dx i dy 



agrees with the traditional ^ when applied to complex-differentiable func- 
tions, yet offers the advantage of applying to the broader class of complex- 
valued real-differentiable functions. Any linear combination of S- and i^- 

J ax 1 ay 
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would work, but the coefficients (|, |) yield the desirable result that 

To take derivatives with respect to z, define 

8 1/8 18 



dz ' 2 \8x i dy J 

and note that |§ = 1 and |§ = 0. For both and J| all the usual differen- 
tiation rules (product rule, quotient rule, power rule) remain valid. 



6 The Laplace Operator 

Parameterize C 2 ~ R A by two complex variables a = x + iy and (3 = z + iw. 
The form of the mixed partials 

d 2 d 2 1 / d 2 d 2 \ 

+ ^ (18) 



da da da da 4 \dx 2 dy 2 
leads us immediately to the complex expression for the Laplace operator 

d 2 d 2 d 2 d 2 



dx 2 dy 2 dz 2 dw 2 
( d 2 d 2 \ 

= 44+41 (19) 



da da 8(3 8(3 



7 Sibling Modes 



Each polynomial of homogeneous degree k in a and (3 alone (no a or (3) 
generates a collection of k + 1 sibling modes. For example, the polynomial 
a 3 — [3 3 generates the four sibling modes 

a 3 -(3 3 
a 2 (3 + (3 2 a 
ap 2 - (3a 2 

p 3 + a 3 (20) 
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Siblings are distinguished by differing twist. Look what happens as we 
trace the value of the original polynomial a 3 — (3 3 (the "first sibling" ) along 
an arbitrary but fixed fiber (at, /3) = (e u ao, e lt f3o). The polynomial evaluates 
to e 3lt (aQ — 0q) on the fiber, so as t runs from to In the polynomial's 
modulus remains constant while its phase runs from to 6n. In other words, 
the first sibling twists three times as we run along a fiber. By contrast, the 
second sibling a 2 /3 + (3 2 a = e^ 2 ' 1 ^ (al(3 + (3 2 a ) twists only once along a 
fiber; the regular variables a 2 and (3 2 contribute two positive twists while the 
conjugated variables j3 and a contribute a negative twist. Similarly the third 
sibling at/3 2 — fta 2 = e^ 1 ' 2 ^ 1 (ao/3 2 — /^a 2 ,) twists minus once along a fiber and 
the last sibling a 3 + (3 3 = e~ 3tt (al + (jq) twists minus three times. 

Kramer [201 generates siblings using what are commonly called "raising 
and lowering operators". Here we will define the same thing but call them 
the negative and positive twist operators 



to emphasize their geometrical significance, which will come in handy in 
Section |H| Algebraically the twist operators' effect is clear: twist _ removes 
a regular variable (via the partial derivatives) and replaces it with a conju- 
gated variable (via the multiplications) for a net decrease of two twists, while 
twist + does the opposite. Geometrically the factorization of a twist operator 
as a dot product (—j3, at) ■ (^, ^) allows an interpretation as a directional 
derivative orthogonal to the fiber. In any case, the reader may verify that the 
positive and negative twist operators take siblings up and down the list (J2(J|) . 
modulo normalization. 

Proposition 01. Modes with different twists are orthogonal. 

Proof. Let / and g be functions with twist m and n, respectively. When 
evaluating the integral (/, g) = L 3 / • g, note that along each fiber (a, (3) = 
(e lt ao, e lt fio) the integral restricts to 



twist. 



twist- 



+ 




(21) 




(22) 
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which equals zero whenever m 7^ n. ■ 

Proposition[7l2. If a symmetry group F C SU(2) preserves a polynomial f , 
then it preserves all the siblings of f as well. 

Proof. Each symmetry 7 = ( z w ~f ) G T defines a rotation of function 
space via the rule [TLy(f)](a, (3) = f[^~ l {a, (3)] for all functions / on S 3 C C 2 
and all points (a, 0) G S 3 . By assumption Ry(f) = f. A quick calcula- 
tion shows that i? 7 commutes with the positive and negative twist operators 
twist± (alternatively, interpreting the twist operators as directional deriva- 
tives orthogonal to the fiber provides deeper geometrical insights, but for 
sake of brevity we will not pursue that interpretation). Thus 

i? 7 [twist ± (/)] = twist± [Ry(f)] = twist± (/) . (23) 

In other words, 7 preserves twist±(f) and therefore preserves all siblings of 
/• ■ 

8 Twist and T-Invariance 

The twist concept sheds some light on why symmetric polynomials turn out 
to be invariant under symmetries of their roots in some cases (Example HJl) 
but not others (Example EJ2). Geometrical considerations often reveal at a 
glance that a given mode cannot possibly be invariant, eliminating the need 
to test the invariance explicitly as we did in Section 

Example [HI 1 . The octahedral group O. An order-2 rotation about the mid- 
point of any of the octahedron's edges lifts to an order-4 corkscrew motion of 
S" 3 (Section^. That order-4 corkscrew motion fixes (setwise) the fiber in S* 3 
corresponding to the fixed edge midpoint on S 2 . Therefore any 0*-invariant 
function must have period 4 along the fiber, meaning its twist must be a 
multiple of 4. Thus the degree-6 polynomial constructed from the octahe- 
dron's six vertices, which has twist 6 along every fiber, cannot possibly be 
0*-invariant. 

Example [HJ2. The tetrahedral group T. As in the preceding example, an 
order-2 rotation about the midpoint of any edge lifts to an order-4 corkscrew 
motion of S 3 , implying that any T*-invariant function must have period 4 
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along the fibers corresponding to fixed points. This would seem to exclude the 
degree-6 polynomial constructed from the tetrahedron's six edge midpoints, 
but in fact that degree-6 polynomial - by construction! - is identically zero 
along the fixed fibers. Thus the degree-6 polynomial could well be (and 
indeed is) invariant under the period-4 corkscrew motion, and in fact turns 
out to be invariant under all of T*. 

On the other hand, the degree-4 polynomial constructed from the tetra- 
hedron's face centers cannot possibly be T*-invariant, because along a fiber 
it's not periodic of order 6, as the order-3 rotation about any vertex would 
require, and is nonzero along half the fixed fibers. 

Example O 3. The icosahedral group I. Section [TO] will construct the de- 
gree 12 polynomial I\2 corresponding to the icosahedron's 12 vertices, the 
degree 30 polynomial /30 corresponding to the 30 edge midpoints, and the 
degree 20 polynomial J20 corresponding to the 20 faces (Equation (J2ZI))- As 
it turns out, luck has smiled on the icosahedral group. 

• The 5-fold rotation about each vertex requires that an invariant poly- 
nomial have period 10 along each fiber corresponding to a fixed point. 
The polynomials J20 and /30 satisfy this criterion automatically because 
they have twist 20 and 30, respectively. While the twist of J 12 is not a 
multiple of 10, its zeros fall at the fixed points of the rotations, so it is 
identically zero along the fibers in question. 

• The 3-fold rotation about each face center requires period-6 invariance 
along the fixed fibers. The polynomials J12 and J30, whose twists are 
multiples of 6, satisfy this automatically. For J 20 the degree is wrong, 
but luckily it's identically zero along the fixed fibers, so it is safe too. 

• The 2-fold rotation about each edge midpoint requires 4-fold periodicity 
along the fiber, which Ji 2 and I20 satisfy by virtue of their twist, and 
/30 satisfies because it's identically zero on the fixed fibers. 

Thus we find no geometrical obstructions to the 7*-invariance of /12, I20 and 
/30 and the method of Section |3] confirms all three to be fully J*-invariant. 
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9 Harmonicity of Modes 

By applying the complex form of the Laplace operator (|19|) to the sample 
siblings (|20J1. the reader may quickly verify that all are harmonic. The fol- 
lowing propositions show that this will always be the case. 

Proposition ®1. Every polynomial in a and (3 alone (no a or (3) is har- 
monic. 

Proof. If / is a polynomial in a and (3 alone, then |^ = = and thus 

Proposition |Hl2. If a polynomial is harmonic then so are all its siblings. 

Proof. The operators V 2 and twist- commute, so if V 2 / = 0, then \7 2 (twist^ f) = 
twist-(V 2 f) = as well, and similarly for twist + . ■ 

Every harmonic function on C 2 ~ M 4 restricts to an eigenmode of the Laplace 
operator V| 3 on S 3 , because 

V|a = r 2 V 2 4 - (r • V)(r • V) - 2(r • V) (24) 

which for a polynomial Q of homogeneous degree k simplifies to 

V§ 3 Q = r 2 V^ 4 g - k(k + 2)Q (25) 

and when Q is harmonic (V^Q = 0) simplifies further to 

V| 3 Q = -k(k + 2)Q. (26) 

10 Modes of the Poincare Dodecahedral Space 

S 3 /P 

For each wavenumber k we will define zero or more base modes. These base 
modes, together with their siblings (Section |7J), will form a complete eigen- 
basis for the Poincare dodecahedral space S 3 /P for the given k. 

For k = the base mode is the constant polynomial 1, which of course 
has no siblings. 
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For k = 12, 30 or 20, construct symmetric polynomials as in Section 
taking P to be the set of vertices, the set of edge midpoints or the set of face 
centers, respectively, of a regular icosahedron, yielding the base modes 



J 12 = a 11 /? + \\a b (3 b - a/3 11 

J 30 = a 30 + 522a 25 /? 5 - 10005a 20 /? 10 - 10005a 10 /? 20 - 522a 5 /? 25 + /? 30 
J 20 = a 20 - 228a 15 /? 5 + 494a 10 /? 10 + 228a 5 /? 15 + /? 20 . (27) 

Still using the method of Section HJ one easily verifies that these three poly- 
nomials are fully invariant under the binary icosahedral group I* C SU(2). 
These three base modes, together with their siblings, form complete eigen- 
bases of dimension 13, 31 and 21, respectively. Historical note: Klein's 
monograph JH] presents these polynomials and notes their invariance under 
the "group of substitutions", but never considers the quotient space S 3 /I*. 
Clifford-Klein space forms weren't conceived until a few years later. 

For wavenumbers k ^ 0, 12, 20, 30, form the base modes as products of 
A2, ho and I 30 . For example, for k = 32 there is exactly one base mode, 
namely I^ho, which together with its siblings forms a basis of dimension 
33. For k = 34 there are no base modes, because 34 cannot be expressed as 
a sum of 12's, 20's and 30's. All products of Iu, ho and /30 are harmonic 
by Proposition Ell and clearly all are invariant under the binary icosahedral 
group I*. 

The wavenumber k = 60 is especially interesting because 60 may be 
written as a sum of 12's, 20's and 30's in more than one way, namely 
60 = 5 • 12 = 3 • 20 = 2 • 30. This gives three candidates for the base 
mode, namely If 2 , -ffo and J| . However, a simple computation reveals a 
linear dependence among the three modes, so any two will span the same 
2-dimensional function space. For sake of concreteness, let us choose 7f 2 and 
i|o, ignoring /| . (In cosmological applications one typically orthogonalizes 
this basis, but for present purposes this isn't necessary.) The two base modes 
l\ 2 and J| , together with their siblings, form a complete eigenbasis of dimen- 
sion 61 + 61 = 122. 

Definition II Oi l. The modulus function mod(m, n) gives the remainder of 
m upon integer division by n. For example, mod(17, 5) = 2. 

Definition 1101 2. The floor function \_x\ gives the integer part of its argu- 
ment x. For example, J = 3. 
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The Poincare dodecahedral space has no modes for odd k, because odd 
modes fail to be invariant under the antipodal map —id El*. To enumerate 
the base modes for arbitrary even k, it is convenient to introduce the notation 
k = k/2. Rather than asking how k may be written as a sum a-12+6-20+c-30, 
we'll ask how k may be written as a sum 



k 



6 + b- 10 + c- 15 

(2 • 3) + b- (2-5) + c- (3-5). 



(28) 



From (|2*H|) it is clear that a = k (mod 5), so in 
: k — 5|_fJ, with notation as in Definitions fTUll 



First consider k modulo 5. 
particular a > mod(fc, 5) = 

andCUl2. Similarly b > mod(£;, 3) = k - 3 [|j and c > mod(£;, 2) = jfe - 2[§J . 
Subtracting off the minimum number of 6's, 10's and 15's leaves 



30 



The quantity 







k 

5 






k - 


- 5 


)• 




k 

5 


+ 


k 
3 


+ 


k 
2 



10 



fe - 2 



15 



k 




k 








+ 




+ 




5 


3 


2 



-jfe 



(29) 



(30) 



tells how many groups of 30 remain after we have subtracted off the minimum 
6's, 10's and 15's. 

If k = —1, then no decomposition k = a ■ 6 + b ■ 10 + c ■ 15 is possible 
(with nonnegative integer coefficients) and the Poincare dodecahedral space 
has no modes of the given order k. If k = 0, then the decomposition is 
unique, consisting of exactly the required minimum number of 6's, 10's and 
15's, and the mode I^I^m along with its siblings spans an eigenspace of 
dimension k + 1. If k > then there will be k factors of degree 60, which 
may be apportioned between If 2 and J| in k + 1 different ways (ignoring J. 

f 2 and /f ). 



because it is a linear combination of i?, and 



2 

30 



Example [T013. When k = 180, then jfe = 90 and k = 3. The k + 1 = 4 base 
modes 

(/i 5 2 ) 3 , (4) 2 (4)), (^X^o) 2 , (4>) 3 , (si) 
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each with its k + 1 = 181 siblings, span an eigenspace of dimension 4 • 181 = 
724. 

Example UM±. When k = 234, then k = 117 = 3- 30 + 2- 6 + 1- 15 and 
there are still the same k = 3 factors of degree 60 as in Example ITU1 1 . but 
now they are supplemented by required minimum factors of degree 12 and 
30. Specifically, I 12 must appear with exponent a = mod(fc, 5) = 2 and 7 3 o 
must appear with exponent c = mod(fc, 2) = 1. No extra factor of I 2 o is 
required because b = mod(/c, 3) = 0. Thus for k = 234 the base modes are 

(/iV3o)(/f 2 ) 3 , (I 2 12 ho)(l! 2 ) 2 (l!o), (/i 2 2 / 3 o)(/ 1 5 2 )(/ 2 3 o) 2 , (AV 3 o)(/ 2 3 o) 3 . (32) 

Together with their siblings they span an eigenspace of dimension 4 • 235 = 
940. 



Theorem 1101 5. For each even k, if we let i range from to 



k/2 



+ 



k/2 



+ 



k/2 



k 
2 : 



(33) 



then the k + 1 base polynomials 



x (a 



20 



(a 11 ^ + lla 6 p 6 - a ^ll)5i+n»d(t,5) 
228a 15 /3 5 + 494a w (3 w + 228a 5 /3 15 + (3 



320N3(fc-i)+mod(f ,3) 



X (a 



30 



+ 522a 25 /3 5 - 10005a 20 /? 10 - 10005a 10 /? 20 - 522a 5 /? 25 + /^mod^,^ 

(34) 

each with its k + 1 siblings (defined in Section^), form a basis of dimension 
{k + l)(k + 1) for the eigenmodes of the Poincare dodecahedral space S 3 /I* . 
There are no modes for odd k. 

Proof. To see that the k + 1 base modes are linearly independent, note that 
they all have different leading terms in the variable a, because a occurs 
with maximal degree 55 in If 2 but maximal degree 60 in J| (Equation (|27Jl ). 
Allowing for the base modes' siblings, the same reasoning applies, but in- 
stead of measuring powers of a alone, measure combined powers of a and 
(3, a quantity which the twist operator (|21j) preserves. This proves that lin- 
ear dependencies cannot occur among the siblings of different base modes. 
Similar reasoning applies within each family of siblings, where differing twist 
expresses itself algebraically as differing powers of the regular variables a and 
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(3 relative to the conjugated variables a and (3. Therefore the full collection 
of base modes and siblings is linearly independent. Technical note: This line 
of reasoning relies on the unique expression of a polynomial in terms of a, a, 
j3 and /3, which follows easily from its uniqueness in terms of the underlying 
real variables x, y, z and w. 

The k + 1 base modes, each with k + 1 siblings (including itself), span a 
space of dimension 



(k + l){k + l) 



( 


k/2 


+ 


k/2 


+ 


k/2 




5 


3 


2 



(k+i) 



(35) 



This agrees with Ikeda's formula for the dimension of the eigenspace |21[ 
Thm. 4.6], proving that our modes span the full eigenspace of S 3 /I*. Note: 
Ikeda computes the spectrum using a dimension counting argument, without 
constructing the modes, so no further comparison is possible. 

As we saw earlier, there are no modes for odd k because such modes are 
not invariant under the antipodal map —id G I*. M 



Reasoning similar to that in the preceding proof shows that for a given 
base mode, sibling i must be the complex conjugate of sibling k — i, for 
% = ... k, up to normalization. In practice the normalization is chosen so 
that corresponding siblings are exactly conjugate to each other, to facilitate 
easy extraction of the real- valued modes. The middle sibling (i = ~) is of 
course its own conjugate up to a factor of ±1. 



11 Modes of the Binary Octahedral Space 

S 3 /0* 

Let us apply the method of previous sections to write down the eigenmodes 
of the binary octahedral space S 3 /0*. Consider in turn the octahedron's 
vertices, its edge midpoints and its faces centers, which according to the 
method of Section |U yield the tentative base modes 

0' 6 = a 5 {3-af3 5 

0' 12 = a 12 - 33a 8 /3 4 - 33a 4 /3 8 + /? 12 

0' 8 = a 8 + 14a 4 /3 4 + /3 8 . (36) 
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The catch is that these base modes are not all 0*-invariant, as Example IHU 
illustrated. More precisely, a quarter turn rotation about any vertex of the 
octahedron lifts to an order 8 Clifford translation of S 3 that preserves 8 
while sending 0' % and 0[ 2 to their negatives. Thus neither 0' 6 nor 0' 12 defines 
a mode of S 3 /0* on its own, but their product 0' 6 0[ 2 and their squares Og 2 
and 0'i2 do. We may safely ignore 0' 2 2 because the mode it represents can 
be recovered later as a linear combination 0' 2 2 
as our base modes 



Of - 108 6 4 . Thus we take 



8 

012 
018 



08 
06 2 

0' 6 0' 12 



a 8 + 14a 4 /3 4 + /J 



a^p 2 - 2a 6 /3 6 + a 2 (3 



a 



10 / 
17 (3 



10 



34a 13 /? 5 + 34a 5 /3 13 - a/3 



?17 



(37) 



Now proceed exactly as in Section [TUl but using {0 8 , Oi 2 , Oi 8 } instead 
of {/i2, ^20) -^3o}- That is, for each even wavenumber k, ask how many base 
modes (if any) may be constructed as products of the {0 8 ,Oi 2 ,Oi 8 }. Still 
imitating Sectional define k = k/2 and ask how k may be written as a sum 



k 



A + b-Q + c-9. 



(38) 



One sees immediately that a = k (mod 3), c = k (mod 2), and 2b = 
k — c (mod 4). The last equivalence reduces to b = ^ (mod 2). Furthermore 

may be 



one may insist that c be or 1, because any occurrence of 0\ 8 



replaced by 2 8 = 0' 2 0' 2 2 = 0' 2 (0' 8 3 - 108 O'i) = Oi 2 (0 3 - 108 O 2 



12/ 



c = mod(A;, 2) transforms the restriction on b to b = mod ( fc - 4 )^ mod ( fc ' 2 ) 
Subtracting from k the minimum multiples of 4, 6 and 9 leaves 



k-3 



12 



HilHHfi: 



+ 



+ 



- k 



k — 2 



Setting 
(mod 2). 



(39) 



and so the quantity 



k 




k 




k 




+ 




+ 




2 


3 


4 



(40) 



tells how many groups of 12 remain. Each group of 12 may be realized as 
either Of or 0\ 2 . 



18 



Theorem II 11 1. For each even k, if we let i range from to 



k/2 



+ 



k/2 



+ 



k/2 



k 
2' 



then the k + 1 base polynomials 



(a S + Ua 4 (3 4 + /?8)3i+mod(!,3) 
X (a W p 2 - 2a 6 p 6 + Q ,2 /3 10)2(fc-i)+(mod(|,4)-mod(| > 2))/2 

x (a 17 /? - 34a 13 /3 5 + 34a 5 /5 13 - ap 17 ) mod ^ > 2 \ 



(41) 



(42) 



eac/i with its k + 1 siblings (defined in Section^), form a basis of dimension 
(k + l)(k + 1) /or i/ie eigenmodes of S 3 /0*. There are no modes for odd k. 

Proof. Same as the proof of Theorem E3 5. ■ 



12 Modes of the Binary Tetrahedral Space 

S 3 /T* 

Let us apply the methods of Sections HO and [TT] to write down the eigen- 
modes of the binary tetrahedral space S 3 /T*. The tetrahedron's vertices, 
edge midpoints and face centers provide tentative base modes 

T' Aa = a 4 + 2ta 2 (3 2 + (3 4 
Tg = a 5 (3-a(3 5 

T ib = a 4 -2ta 2 (3 2 + (3 4 . (43) 

A one-third rotation about any vertex of the tetrahedron lifts to an order 6 
Clifford translation of S 3 that preserves Tg but sends T' Aa \— > £ T' Aa and T' 4b i— > 
£ T' Ab , where £ = | + is a third root of unity Therefore the group T* 
preserves the modes 

T 6 = Tg = a 5 f3 -af3 5 

T 8 = T ia T Ab = a 8 + Ua 4 p 4 + P 8 (44) 

T 12 = (Ti a 3 + Ti 6 3 )/2 = « 12 - 33« 8 /? 4 - 33« 4 /? 8 + /? 12 . 

We may ignore the complementary order 12 mode (T^ 3 — T' ib )/2 because it 
is a multiple of T 6 2 . Note that the tetrahedral modes (jT4j) coincide exactly 
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with the tentative octahedral modes ([SB]), reflecting the close ties between 
O* and its index 2 subgroup T*. 

Let us construct the modes of S 3 /T* as in the preceding sections, by 
taking products of {T 6 ,T 8 ,Ti 2 } to establish the base modes and then taking 
all their siblings. With k = k/2 as before, ask how k may be written as a 
sum 

fc = a- 3 + 6- 4 + c-6. (45) 

One sees immediately that a = k (mod 2) and b = k (mod 3). Subtracting 
from k the minimum multiples of 3 and 4 leaves 



k - 



k-2 
6 



k 
2 

+ 



3 - ( jfe - 3 
+ 



(46) 



and so the quantity 



k 




k 




k 




+ 




+ 




2 




3 


3 



(47) 



tells how many groups of 6 remain. Each group of 6 may be realized as either 
T 6 2 or T 12 . We may safely restrict T 8 to the minimum b = 0, 1 or 2 occur- 
rences (according to the value of k mod 3), because additional occurrences 
may be swapped out via the dependency T| = Tf 2 + 108Tg. 



Theorem I12L 1. For each even k, if we let i range from to 



k/2 



+ 



k/2 



+ 



k/2 



k 
2' 



(48) 



then the k + 1 base polynomials 



{a 5 fJ - c^5)2i+mod(f ,2) ( a 8 + ^4^4 + ^8) mod( | ,3) 



x (a 



12 



33a 8 /3 4 - 33a 4 /3 8 + /? 12 ) fc " 



(49) 



eac/i its k + 1 siblings (defined in Section^, form a basis of dimension 
(k + l)(k + 1) for the eigenmodes of S 3 /T* . There are no modes for odd k. 

Proof. Same as the proof of Theorems 5 and II 11 1. ■ 
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13 Modes of the Binary Dihedral Spaces S 3 /D* 

Still following the method of the preceding three sections, let us write down 
the eigenmodes of the binary dihedral space S 3 /D^ n , where D* m is the binary 
dihedral group of order 4m. The dihedron's vertices (sitting at m th roots 
of 1), edge midpoints (at m th roots of —1) and face centers (at {0, oo}) 
provide the tentative base modes 

D' v = a m -(3 m 
D' e = a m + (3 m 

D' f = a(3. (50) 

A half turn about the vertex at 1 lifts to the order 4 Clifford translation 
(a, /3) i-> = (i(3,ia), which sends D' v ^ -i m D' v , D' e ^ +i m D' e , and 

D'j i— »• — D'j. Similarly a half turn about the edge midpoint at e™ lifts to the 
order 4 Clifford translation (a,f3) i— > ( °_ i2L te m ) (%) = (ie™(3, ie~a), which 

\ie rn / 

sends D' v i-> +i m D' v , D' e ^ -i m D' el and D' f ^ -D'f. These two half turns 
generate the full group. 

Keeping in mind that the product D' v D' e = a 2m — f3 2m is invariant when 
m is odd, while (D' e ) 2 — 2(D'f) m = a 2m + (3 2m is invariant when m is even, 
we may combine the tentative base modes (|5U|) into the final base modes 

D A = a 2 (3 2 

D 2m = a 2m ±P 2m (51) 
D 2 m+2 = a(3(a 2m TP 2m ) 

with the upper sign choices for m even and the lower sign choices for m odd. 
In both cases the modes are fully invariant under the two half-turns and 
therefore under all of D^. 

For each wavenumber k ask how many independent base modes may be 
constructed as a product 

DtD 2m D2 m+2 . (52) 

We may safely insist that c be or 1, because the square D 2m+2 may be 
replaced by D\ m+2 = D A (D 2 m =p 4DJ 1 ). To facilitate an inductive approach, 
let Bk denote the set of modes constructed for wavenumber k. For k = 
the only mode is the polynomial 1, so Bq = {1}. For k = 2 there are no 
base modes when m > 1, so B 2 — {} is empty, but there is a base mode 
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B 2 = {a 2 — /3 2 } when m = 1. Now proceed inductively. Assuming we've 
already constructed the set -Bfc-4, let B k consist of the polynomials in B^-^ 
each multiplied by an additional factor of D4, along with an additional mode 

^ k \ ^ fc-(2m+2) s 

D^m ^ 2m divides k or D 2m 2m D2 m +2 if 2m divides k — 2 (or both in 
the special case m = 1). To prove that the polynomials in Bk are linearly 
independent, note that each one has a different highest power of a. 

A few concrete examples illustrate the iterative — and ultimately very 
simple - nature of the construction. The reader is encouraged to carry out the 
construction, as defined in the preceding paragraph, for the cases m = 1,2,3 
and compare his or her results to those given in the following three examples. 

Example I13L 1. In the seemingly (but not really) exceptional case m = 1 
the construction proceeds as follows: 

Bo = {1} 

B 2 = {a 2 -f3 2 } 

B 4 = {a 2 p 2 , (a 2 -p 2 ) 2 , a(3(a 2 + (3 2 )} 

B 6 = {a 2 (3 2 (a 2 -(3 2 ), (a 2 -(3 2 ) 3 , (a 2 -(3 2 )af3(a 2 + (3 2 )} (53) 
B 8 = {a 4 (3\ a 2 p 2 (a 2 - (3 2 ) 2 , a 3 (3 3 (a 2 + (3 2 ), 
(a 2 -(3 2 )\ (a 2 - (3 2 ) 2 a(3(a 2 + (3 2 )} 



Example I13L 2. For m = 2 the construction yields 

Bo = {1} 

B2 = {} 

B A = {a 2 (3 2 , a 4 + (3 4 } 

B 6 = {af3(a A - /? 4 )} 

B 8 = K/3 4 , a 2 p 2 (a 4 + p 4 ), (a 4 + (3 4 ) 2 } 

B w = {a 3 (3 3 (a 4 - /? 4 ), (a 4 + f3 4 )a[3(a 4 - (3 4 )} 

B 12 = {a 6 p 6 , a 4 p 4 (a 4 + p 4 ), a 2 p 2 (a 4 + p 4 ) 2 , (a 4 + p 4 ) 3 } 
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Example I13L 3. For m = 3 the construction yields 

{1} 
{} 

{a 6 - (3 6 } 

{a 2 {3 2 (a 6 -(3 6 )} 
{a 6 /? 6 , a 3 (3 3 (a 6 + {3 6 ), 



Bo 
B 2 
B 4 
B 6 
B 8 
Bio 

B\2 



(55) 



[a 



For all k > and all m > (including m = 
shows the number of base modes WA to be 



1) a straightforward induction 



k/2 



+ 



fc/2 



+ 



k/2 



m 



k - + i. 

2 



(56) 



This formula looks superficially different from the result of Ikeda's dimension 
counting argument [21] , which states separate formulas for k = (mod 4) and 
k = 2 (mod 4), but turns out to be equivalent. Note, however, that Ikeda's 
Theorem 4.3 needlessly excludes the case m = 1 and accidently excludes the 
eigenmodes for k = 2m when m is odd. 

The space S 3 /D^ n admits no modes of odd wavenumber k because such 
modes fail to be invariant under the antipodal map —id e D* m . 

Theorem I13L 4. For each m > and each even wavenumber k, the polyno- 
mials 

(a 2 f3 2 ) a (a 2m + (3 2m ) b {a(3{a 2m - (3 2m )) c if m is even 
or (57) 
{a 2 f3 2 ) a {a 2m - f3 2m ) b {a(3{a 2m + f3 2m )f if m is odd 

satisfying a > 0, b > 0, c G {0, 1} and 4a + 2mb + (2m + 2)c = k, together 
with all their siblings (defined in Section^, form a basis for the eigenmodes 
of S 3 /D* m . The basis has dimension 

-± + l\(k + l). (58) 



( 


k/2 


+ 


k/2 


+ 


k/2 




2 


2 


m 



and the polynomials may be easily constructed by the iterative procedure de- 
scribed above. There are no modes for odd k. 



Proof. Contained in preceding discussion. 
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14 Modes of the Lens Spaces L(p, 1) 



The "binary cyclic groups" Z* give only even order lens spaces, so rather 
than trying to lift group elements and modes from S 2 to S 3 as we did in 
Sections ITUlTEfl let us instead work in S 3 C C 2 directly The lens space 
L(p, 1) is generated by the Clifford translation 




which immediately tells us that the group preserves the mode ot a (3 b if and 
only if a — b = (mod p). Thus for a given wavenumber k (odd or even) the 
number of base modes is simply the number of ways to choose nonnegative a 
and b satisfying a + b = k and a — b = (mod p). Examining these constraints 
on the integer lattice in the (a, b) plane gives the dimension counts stated 
in Theorem [TH 1 . These counts agree with previously obtained real- variable 
eigenmodes for lens spaces [22J. 



Theorem 1141 1. For each p > and k > (odd or even), the monomials 

a a fJ b (60) 

for a>0,b>0,a + b = k and a — b = (mod p), together with all their 
siblings (defined in Section^, form a basis for the eigenmodes of L(p, 1) = 
S 3 /Z p with dimension 

{^lp\ + lj (&+ 1) if p is even and k is even 
if p is even and k is odd 

(2L4J + !)(*+!) if pis odd and. is even <«' 

2 (L^J + l) (Jfe + 1) if P is odd and k is odd 

In the special case of the "binary cyclic groups" S 3 /Z* = L(2n, 1) the dimen- 
sion equals 

+ 1)- (62) 



( 


k/2 


+ 


k/2 






n 




n 





for even k and zero for odd k, in analogy with Theorems \1(A5. [771 1.M2L1 and 

mi 



Proof. Contained in preceding discussion. 
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15 Modes of the 3-Sphere S 3 

For sake of completeness, let us note that the 3-sphere fits comfortably within 
the framework of the preceding sections. 

Theorem I15L 1. For each k > (odd or even), the monomials 

a a (3 b (63) 

for a > 0, b > and a + b = k, together with all their siblings (defined in 
Section^, form a basis for the eigenmodes of S 3 with dimension (k + l) 2 . 

Proof. The polynomials a a f3 b and their siblings are linearly independent, are 
harmonic by Propositions 01 and 02, and meet the well-known dimension 
(k + l) 2 of the mode space of S 3 . ■ 



16 Future Work 

The present article provides explicit polynomials for the eigenmodes of all 
globally homogeneous spherical 3-manifolds. Four directions present them- 
selves for future work. 



16.1 Find a unified approach 

Mathematically, the similarity of the base mode dimension formulas (|33|). 
(|4ip . (|48p. (Jo"S|) and ()62j). along with their obvious dependence on the under- 
lying pqr triangle group, is too striking to ignore. While this could be merely 
an effect of dimension counts, one suspects a single unified proof might be 
possible. Here is a plausible approach, for all cases except the odd order lens 
spaces. Starting with one of the base modes, whose twist equals its (even) 
degree k, apply the twist- operator k/2 times to obtain a sibling of twist 0. 
Because the twist is zero, the sibling projects down from a mode of S 3 to 
a well-defined mode of S 2 . Elementary considerations imply that the mode 
is pure real or pure imaginary, so after possibly multiplying by % we may 
assume it is pure real. Even more elementary considerations imply the result 
is harmonic on S 2 . The idea for a unified proof is to reverse-engineer this 
process. First construct polynomials on S 2 forming a basis for the modes 
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invariant under a pqr triangle group, then lift those modes to twist modes 
of S 3 , and finally apply both the twist- and the twist + operator k/2 times 
each to obtain all siblings. 

16.2 Orthonormalize the basis 

The theorems proved in this paper provide a basis for the space of eigenmodes 
of every homogeneous spherical 3-manifold. CMB simulations, however, re- 
quire not just an arbitrary basis, but rather an orthonormal basis. The author 
has tried to find a simple and elegant way to normalize the given bases, but 
so far without success. The problem seems to require a fundamentally new 
insight. As a fallback plan, numerical methods would work, but with a heavy 
cost in terms of slower computation and diminished accuracy. 

16.3 Transfer to spherical coordinates 

The modes found in this paper are wonderfully simple: they are explicit 
polynomials of fairly low degree in only two variables. For CMB simulations, 
however, we need to relate these modes to the last scattering surface, which in 
effect requires transferring the modes to spherical coordinates. The change- 
of-coordinates must be done quickly and accurately; otherwise the explicit 
polynomials will not be significantly faster or better than numerical methods 
(e.g. using more than 10000 eigenfunctions [H]). 

16.4 Improve CMB simulations 

Cosmologically, the next step is to apply the polynomials we have found to 
streamline the cosmological simulations. The late Jesper Gundermann be- 
gan this process and reported improved accuracy, the elimination of small 
extraneous imaginary quantities from the calculation, and the possibility of 
extending the cosmological simulations to higher wavenumbers k. A full ap- 
plication of the polynomials will require satisfactory solutions to the problems 
discussed in Sections 116.21 and 116.31 
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